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We analyze, via Imry-Ma scaling arguments, the strong disorder phases that exist in low dimen- 
sions at all temperatures for directed polymers and interfaces in random media. For the uncorrelated 
Gaussian disorder, we obtain that the optimal strategy for the polymer in dimension 1 + d with 
< d < 2 involves at the same time (i) a confinement in a favorable tube of radius Rs ~ L"^ with 
Us = 1/(4 - d) < 1/2 (ii) a superdiffusive behavior R L" with i/ = (3 - d)/(4 - d) > 1/2 for the 
wandering of the best favorable tube available. The corresponding free-energy then scales as _F ~ 
with = 21^ — 1 and the left tail of the probability distribution involves a stretched exponential of 
exponent 77 = (4 — d)/2. These results generalize the well known exact exponents v = 2/3, cj = 1/3 
and rj = 3/2 in d = 1, where the subleading transverse length Rs ^ L}^^ is known as the typical 
distance between two replicas in the Bethe Ansatz wave function. We then extend our approach 
to correlated disorder in transverse directions with exponent a and/or to manifolds in dimension 
D + d = dt with < D < 2. The strategy of being both confined and superdiffusive is still optimal 
for decaying correlations (a < 0), whereas it is not for growing correlations (a > 0). In particular, 
for an interface of dimension (dt — 1) in a space of total dimension 5/3 < dt < 3 with random-bond 
disorder, our approach yields the confinement exponent vs ~ {dt — l){3 — dt) / (5dt — 7). Finally, we 
study the exponents in the presence of an algebraic tail l/V^"'"^ in the disorder distribution, and 
obtain various regimes in the (/i, d) plane. 

I. INTRODUCTION 

Directed polymers in random media have attracted a lot of interest for many years, either as interesting disordered 
models or in relation with stochastic growth models . For polymers in dimension 1 -I- d described in the continuum 
limit by the partition function 

' Pf(s)e- -lo HST-'^ lo 'isV{s,r{s)) 

with an uncorrelated Gaussian random potential 



Vis,f}Vis',f') = 6{s^s')S'^{r-r') (2) 

the phase diagram is the following : in dimension d > 2, there exists a phase transition between a free phase at 
high temperature and a pinned phase at low temperature : this phase transition has been studied numerically 

in d = 3 Q, exactly on a Cayley tree and on hierarchical lattice On the contrary, in dimension d < 2, there is 
no free phase, i.e. any initial disorder drives the polymer into a strong disorder phase. The marginal dimension d — 2 
has been controversial and deserves a special discussion A strong disorder phase is characterized in particular by 
two exponents u and v for the free energy F and the transverse length scale R 

F{L) - L"^ (3) 
R{L) ^ (4) 

with the expected scaling relation w — 2v—\ 0. In 1 + 1, these exponents uj and v are exactly known to be w = 1/3 and 
^ J/ = 2/3, because in c? = 1, some "miracles" happen in various methods : via the mapping towards a damped Burgers 
■ ■ ■ ' equation with random forcing, there exists an exact steady-state distribution that fixes the values of the exponents "^l; 

within the replica framework, there exists exact Bethe Ansatz solutions, that have been studied either in unbounded 
space 1^ ^1 , or in bounded space ; there exists an exact combinatorial solution at zero-temperature i as well 
as other exact results via the correspondence with stochastic growth models 13] in the same KPZ universality class. 
So there are plenty of reasons why the exponents are exactly 1/3 and 2/3 in d = 1. However, various questions are 
still open or under debate concerning the generalizations of these exponents in various ways, namely : 

(i) in other transverse dimensions d 

(ii) in the presence of transverse correlations 

(iii) for manifolds and interfaces of higher internal dimensions D 

(iv) for various initial disorder distributions, presenting for instance algebraic tails. 
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In this paper, our aim is to present Imry-Ma scaling arguments that allow to analyze these various generalizations 
in a unified framework, by a proper identification of the underlying optimal strategy in each case. The paper is 
organized as follows. 

In Section (nj we recall the two 'local' Imry-Ma arguments proposed in for the directed polymer in 1 + c? 
dimensions in favorable and unfavorable regions. In SectionlTTIl we propose a global optimization mechanism between 
the energy gained by a confinement in a favorable 'tube' and the global elastic energy to find the best favorable tube 
available. This strategy fixes a confinement exponent, a global wandering exponent, as well as a free-energy exponent 
and the form of the left tail for the free energy probability distribution. In Section llVI we generalize our approach 
to other correlations in transverse directions and/or to other internal dimensions : we discuss in particular the cases 
of directed polymers with correlated transverse disorder, of interfaces with random bond disorder and of interfaces 
with random field disorder. In all these cases, we find exponents in agreement with the replica-scaling analysis by 
Zhang |2lj| , and we thus discuss the dictionary between our approach via Imry-Ma argument in one disordered sample 
and the Zhang analysis in replica space. Finally in Appendix ^ we generalize our approach for directed polymers 
in dimension 1 + d in the presence of an algebraic tail 1/^^+^ in the disorder distribution, and we obtain various 
regimes in terms of the two parameters (/i, d). 



II. TWO LOCAL IMRY-MA ARGUMENTS FOR FAVORABLE AND UNFAVORABLE REGIONS 

In this Section, we recall in details the two Imry-Ma arguments proposed in [l^ which constitute the basis of all 
our discussion. (In Il4| , these arguments were given to interpret the two types of solution found via the disorder- 
dependent variational method). As in other contexts, the Imry-Ma argument |l5l | begins with the evaluation of the 
typical energy associated to the disorder in a certain volume. Here, for a polymer of length I and transverse length r, 
the dimensional analysis of the correlator ||2Jl yields the following scaling for the typical random energy 

(j^dsV{s,r{s))\ ^±u\f^ (5) 

\ / typ 

where it is a random variable of order 1. Regions with u — u- > correspond to globally "unfavorable regions", 
whereas regions with u = — u+ < correspond to globally "favorable regions" . 

In other contexts, such as random-field Ising models the Imry-Ma argument then consists in comparing the 
energy cost in creating domain walls with the typical energy gained by taking advantage of the favorable fluctuations 
of the disorder. Here in the polymer context, the energy coming from the disorder has to be compared with the 
entropy cost, that can take two different forms '■ for a swollen polymer r ^ Z^^^, the entropy cost consists in an 
elastic term Tr"^ /I , whereas for a confined polymer r <C V'^'^, the entropic cost consists in a confinement term Tl/r^. 
These t wo p ossibilities lead to two different Imry-Ma arguments, that can be associated to unfavorable and favorable 



;wo_p' 

regions [lj| as we now explain. 

A. Imry-Ma argument with the elastic term for "unfavorable regions" 

The free-energy of an unfavorable region of length I- and transverse length r„ is the sum of the elastic term Tr^^/l-, 
that represents an entropic cost, and the energy cost u-^J^- from the unfavorable fluctuation of the disorder Q 

/-~^7^ + «-W-^ (6) 
The minimization with respect to r_ yields, after dropping numerical prefactors 

— with i._ = -^ (7) 

T ) - 4 + d ^ ' 

in particular v^{d = 1) = 3/5. The corresponding scaling for the free energy © of this unfavorable region reads 

f ^ Tifiiu^r- with = 2j._ - 1 = ^— ^ (8) 
" " 4 + d ^ ' 
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in particular u}-{d = 1) = 1/5. These exponents v- and lo- actually correspond to the direct dimensional analysis of 
the initial Hamiltonian, and are usually called "Imry-Ma exponents" or "Flory exponents" in the more general context 
of interfaces and manifolds in random media [llllllIllllSElEllla. In dimension d ^ 1, these exponents 
= 3/5, CO- = 1/5) are also the exponents predicted for the full polymer by the Replica Gaussian variational Ansatz 
with replica symmetry breaking in contrast with the correct exponents {v — 2/3, uj = 1/3) for the full polymer 
found by the replica symmetric Bethe Ansatz solution 9J. 

Here, we stress that the above Imry-Ma dimensional analysis should not a priori be applied blindly to the full 
polymer, but only to the unfavorable regions. Our conclusion for the moment being is thus the following : if the 
polymer has to cross an unfavorable region ii_ < 0, it will behaves as follows when the dimension d varies : 

(i) for < d < 2, the polymer will adopt a wandering exponent z/_ ~ > 1/2, and the free energy will have for 

exponent w_ = > 0. 

(ii) for the marginal case d = 2, the wandering exponent reaches the free value i/^ld = 2) = 1/2 and the free energy 
exponent vanishes u;-(d = 2) = 0. A more refined analysis thus becomes necessary. 

(iii) for d > 2 the above Imry-Ma argument that would yield v- < 1/2 breaks down, since for a confined polymer, 
the elastic free energy r'^/l- has to be replaced by the confinement free energy l-/rt. However, in this case, the free 
energy is minimum in the limit r_ oo, that does not correspond to a confined configuration. So at the level of this 
scaling analysis, the only consistent possibility for c? > 2 is that the polymer will keep its free exponent j/_ — 1/2 
corresponding to a finite elastic free energy, and the disorder potential then corresponds to a subleading term of order 
/(2-d)/2^ which is what happens in the high temperature phase. 



Imry-Ma argument with the confinement term for "favorable regions" 



The free energy of a favorable region (Z+,r+) is the balance between the confinement term TZ+/r^, representing 
le 



the the entropy loss due to the confinement, and the energy gain from a favorable fiuctuation of the disorder 



/+-ri±-u+./i± (9) 

The minimization with respect to yields, after dropping numerical prefactors, 

— 11+ with 1^+ = (10) 

in particular ^^(d = 1) = 1/3. The corresponding scaling for the total free energy of this favorable region reads 

f+ ^ -T-'^''+u*//+lX'^ with = 1 - 2i/+ = |— ^ (11) 

in particular uj^(d = 1) = 1/3. In contrast with the exponents (j/_,cij_) coming from a direct dimensional analysis 
of the Hamiltonian, the exponents (z^_(-, w-j-) take into account the physical idea that it can be better for the polymer 
to remain confined in a region to take advantage of favorable fiuctuations of the disorder. To our knowledge, these 
exponents (j/+,cij+) have not been considered previously, except in jT3 | where they have been introduced. 

Our conclusion for the moment being is thus the following : if the polymer has to cross a favorable region u+ > 0, 
it will behaves as follows when the dimension d varies : 

(i) for < d < 2, the polymer will adopt a confinement exponent = < 1/2, and the free energy will have 

for exponent w+ — > 0. 

(ii) for the marginal case d = 2, the confinement exponent reaches the free value I'+id — 2) — 1/2 and the free 
energy exponent vanishes a;+(d = 2) — 0. A more refined analysis thus becomes necessary. 

(iii) for d > 2 the above Imry-Ma argument that would yield > 1/2 breaks down, since the polymer is not 
confined anymore. If one replaces the confinement term l+/r'^ by the elastic term r'^/l+, the total free energy will be 
minimum for r+ — > corresponding to a confined configuration. So at the level of this scaling analysis, exactly as in 
unfavorable regions, the only consistent possibility for d > 2 is that the polymer will keep its free exponent — 1/2 
corresponding to a finite elastic free energy, and the disorder potential then corresponds to a subleading term of order 
l{2~d)/2^ which is what happens in the high temperature phase. 
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C. Discussion 



In this Section, wc have described via Imry-Ma arguments what typical scahngs should be expected from a polymer 
that is obliged to go through a given favorable region or through a given unfavorable region. In dimension < d < 2, 
this analysis yields two sets of non-trivial exponents {v_,l^_) and {i'+,uj+) for the two types of regions, whereas 
for d > 2, the only self-consistent exponents in the above Imry-Ma scaling analysis are the exponents of the high- 
temperature phase. This suggests that the pinned phase existing in dimension d > 2 at low temperature is very 
different in nature from the physics in dimension d < 2 and requires a different type of analysis. In the following, we 
will thus only consider the cases < d < 2, where the disorder is strong at all scales and changes the free exponents 
both in favorable and in unfavorable regions. 



III. STRUCTURE OF THE FULL POLYMER OF LENGTH L IN DIMENSION < d < 2 

In this Section, we consider the standard situation of a polymer of length L whose origin is fixed and whose end- 
point is free. We discuss what is the best strategy to obtain a minimum free energy, in terms of the favorable regions 
described in previous Section. 



A. Global optimization on scale L and exponents 



The simplest strategy that seems optimal at large scale L is the following : the polymer will try to find a favorable 
region of length ~ L, of transverse length i?+ ^ and of free energy ^ —L'^+. The only degree of freedom 
available to find this very favorable region is the global orientation Rq ~ , with respect to the horizontal line, of the 
tube of radius starting from the origin forming the favorable region. To find the best favorable region available, 
the polymer can afford a global elastic cost TRq/L that is at most of the same order of magnitude of the free energy 
F+ ~ L"+ of the favorable region it is looking for. The balance between these two terms yields the following global 
transverse distance 



and the corresponding free energy 



Rg with ^ = = ^ (12) 



with to = ujjL = (13) 

4-d 



These two exponents are thus the generalizations in dimension < < 2 of the well known exact exponents v — Iji 
and aj = l/3inc?=l Moreover, our description also yields the subleading transverse length scale 

i?+ - with 1/+ = - (14) 
4-d ^ ' 

representing the radius of the "tube" of the favorable region, that generalizes the transverse length scale L^/^ in- 
troduced in [2&] to characterize the size of a "family", i.e. paths having free energy differences of order 0(1). This 
subleading transverse length scale L^l^ was also interpreted in as the typical distance between two replicas in the 
Bethe Ansatz replica wavefunction 0, whereas the scale L?^^ represent the displacement of all replicas as a whole. 

The physical meaning of the present Imry-Ma scaling analysis is thus the following : the configuration of the 
polymer is determined by a global optimization mechanism at the largest scale : the polymer chooses the best "tube" 
of radius i?+ ^ L'^+ among all tubes available labeled by the global orientation p = Rj^ / L'^ defined by the transverse 
distance i?i of the end point. The number of different tubes available for the choice thus scales as 



Ri 



N^-^^L-' with ^ = d(^_^+)^^_^ (15) 



in particular 7(d = 1) = 1/3. This number is large enough to find a 'good' tube in a arbitrary sample. 
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B. Tail of the distribution of the free energy 



The Imry-Ma analysis for favorable regions can be used to study the tail of the probability distribution of the 
rescaled free energy. Indeed, the asymptotic behavior of the random variable u+ in the Imry-Ma argument (|5I9|I is 
expected to follow to the Gaussian distribution 

P{u+) ~ e-< (16) 

U-^ — >oo 

The same idea has been already used in the context of the random field XY model 0|, and in the context of a 
disordered heteropolymer , where it was shown to be in full agreement with a disorder dependent real-space RG 
analysis (see the Appendix of [l^). In the present context, we stress that the Gaussian tail (|16ll is valid for an initial 
Gaussian disorder, whereas the presence of an algebraic tail of arbitrary order in the initial disorder will generate a 
different tail, as discussed in Appendix 

Here, for a Gaussian initial disorder, the Gaussian tail H16|l yields, via the change of variables the following 
decay for the probability distribution of the rescaled free energy /+ = / L'^ 

n(/+) oc T'^/^(|/+|)^-e--^'^l/+l- with = = i-i^ (17) 

At the level of exponents where we work, since L-|_ ^ L, the tail of the probability distribution P(/) for the rescaled 
free energy f ^ F/ L'^ oi a, polymer of length L will thus be given by the same form (|17|l 

V{f) ^ - V+{f) (18) 

In particular, the exponent in the exponential in one dimension is r]+(d = 1) = 3/2, a value that agrees with the 
replica-scaling predictions 0, |^ and with the numerical simulations [23] ■ 

However, to be fully consistent with the approach we propose, we should take into account that the polymer actually 
chooses the best tube among of large number N (|15|l of independent tubes. In this interpretation, the random variable 
M+ is not just a random variable drawn with a distribution having the Gaussian tail H16|l. but it is the maximal value 
Umax drawn among N independent variables, i.e. its distribution reads 



PN{Umax) = NP{u„ 



+ 00 



duP{u) 



N-l 



N- 



AT —u^ —N r + ^ due ^' /I n\ 



In particular, Umax grows logarithmically in N and thus in L (|15f) 



Umax viniV - VlnJ (20) 
This would lead Hll|) to a logarithmic correction to the exponent for the free-energy (|13|l . i.e. 

F^-L-ut:i.,^-L-i\nL)^ (21) 

Is this logarithmic correction a reality or an 'artefact' of our interpretation? On one hand, the comparison with 
the results of most other studies on the subject suggests that this logarithmic correction is spurious. On the other 
hand, within our approach, it is not clear to us why this logarithmic correction should be disregarded. In particular, 
if the initial disorder distribution is not Gaussian but presents an algebraic tail of index (1 -f /i), it is precisely this 
mechanism of choosing the best variable u^ that opens the possibility of obtaining different global exponents even if 
the variance is finite /i > 2, as discussed in Appendix ^ We note moreover that the presence of some logarithmic 
factors coming from extremal statistics has already been proposed and numerically studied for the directed polymer 
in 1 -|- 1 p3 |. as well as in another context involving polymers in random media |29j |. In conclusion, even if the exact 
solution at zero temperature |l2l | has no logarithmic correction, since we are not aware of exact solutions at non-zero 
temperature, and since the limit of zero-temperature cannot be discussed within our approach (see the discussion 
below in llll C|l . it seems that the presence of this logarithmic correction for the free energy at non-zero temperature 
is a possibility that cannot be completely ruled out, at least to the best of our present knowledge. If there is a proof 
in the future that there is no logarithmic correction, this would probably mean that the 'best' tube is not simply the 
tube having the maximal rescaled variable u, but the tube having the best structure on smaller scales than the global 
scale. 
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C. Remarks on the zero-temperature limit 

In this paper, we have considered that the temperature appears only in front of the random potential in the partition 
function Q and not in front of the Wiener measure for the Brownian paths. In this case, the elastic term is an entropy 
coming from the probability 1^ to be at distance K in time L for a Brownian motion. In particular, this elastic 
term is not present at T = 0, where the problem on the hypercubic lattice in (1 + d) consists in finding the best path 
among the (1 + d)^ possible paths, i.e. on the lattice, a path is either allowed or forbidden, there is no elastic energy 
at T = 0. This corresponds to the usual model for numerical simulations on directed polymers at T = 0. However, 
many authors are also interested by the models where the elastic term jL is an energy, i.e. in the partition function 
(|H) , the temperature appears not only in front of the random potential but as a global factor in front of the two terms 
in the exponential In this context, the problem at T = consists in finding the best path that optimizes the sum 
of the elastic energy and the random energy. 

However, in both cases, within our approach, the limit of zero-temperature T — > is very singular, because the 
entropy due to confinement that played a crucial role at T ^ disappears at T = 0. Let us briefly see what happens 
as T — !■ in the Imry-Ma arguments energy/entropy that we have discussed. For unfavorable regions, the prefactor 
of the transverse scale r_ JTI) diverges and the prefactor of the free-energy © vanishes : indeed, at T = 0, the 
free-energy only contains the random energy from the disorder (and no entropic elastic term anymore) : its 
minimization corresponds to a transverse scale r_ as big as possible, i.e. on a lattice r_ ~ /_ that leads to an energy 
e_ ~ U-tli Pqj, favorable regions, the prefactor of the confinement scale r+ (|10() vanishes, whereas the prefactor 

of the free-energy diverges : indeed, at T = 0, the free-energy (j^J only contains the random energy from the 
disorder (and no entropic confinement cost anymore) : its minimization corresponds to a confinement scale r+ as small 

1/2 

as possible, i.e. on the lattice r_|_ ^ 1 corresponding to a unique path, that leads to an energy of order e+ ~ —uj^Ij^ . 

In conclusion, our description with Imry-Ma arguments cannot be used to understand the zero-temperature limit, 
even if there are very direct relations between T 7^ and T = : in d = 1, the wandering exponent of the best 
tube at finite temperature = 2/3 is the same as wandering exponent of the best path at zero temperature, and the 
subleading transverse scale -R+ ^ L^l^ that represents in our approach the confinement scale at non-zero temperature, 
had been previously identified in zero-temperature numerical simulations as the typical scale for the first excited states 
of finite energy above the ground state Hq, |33 • 

D. What is the sub-structure of the polymer at smaller scales ? 

From the point of view of the Imry-Ma scaling analysis proposed in this paper, it is clear that the global exponents 
Lo and v (|12|l are completely constrained by a global optimization mechanism at the biggest scale. However, once 
the best global tube has been chosen, it seems natural to expect that the polymer has a "cascade" of optimizations 
to make on smaller and smaller scales within the large scale constraints. In particular, we expect that the extensive 
contribution to the free-energy comes from the small scales, since the polymer has to gain a finite contribution at 
each step on average. However, we also expect that the polymer cannot avoid frustration on all scales, i.e. it will be 
obliged sometimes to cross "unfavorable regions" characterized by the exponents and uj- (|7I8|I . A first interesting 
question is : what is the scale L_ of the biggest unfavorable region inside the globally favorable tube? A more general 
question concerns the hierarchical organization on smaller and smaller scales. Indeed, the ultrametric tree structure 
of local optimal paths [s^ [s^ as well as Monte-Carlo simulations at finite temperature Q suggest that they are 
favorable and unfavorable regions on various scales. In conclusion, it would be very interesting to have a statistical 
description of the substructure of the global tube, but this goes beyond the present work. 

E. Discussion of some consequences 

The idea of favorables tubes where the polymer remains confined at finite temperature, with a confinement radius 
much smaller than the wandering scale ^ gives a more precise meaning to the notion of 'states' developed in 
the droplet scaling theory ^31, 32] . It is thus interesting to mention briefly some important consequences in dimension 
< d < 2 with the values of the global exponents that we have obtained. (It would be of course very interesting to 
give a more precise characterization of the 'states' at low temperature in dimension d > 2, but this is left for future 
work) . 
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1. Statistics of the effective random potential for the end-point 
The effective Hamiltonian seen by the free end-point r — r{L) can be decomposed into [sil l32 



It 



i^e//=T- + L"<I>^p.-J (22) 

where the first term represents the elastic free energy and the second term an effective random potential that has 
been rescaled with the global exponents, and whose statistics has to be elucidated. The rescaled effective potential 
<I>(/3), has been exactly determined in d = 1 : it is an 'Airy process' that behaves locally as a random walk ^/p 
as p ^ and that saturates towards a constant at large distances p ^ oo, in agreement with the previous numerical 
study ^(|. More generally, in dimension d, the rescaled effective potential is expected to be independent of L at 
short distance p — *■ |31| , and thus the exponent cr defining the power-law behavior of the effective potential at short 
distances 

<&(/?) cx IpT (23) 
is not a new exponent, but is a function of the two basic exponents |3ll l82| 

(24, 

With the values obtained before H12I13|I . we thus obtain in dimension < < 2 

2 - d 
- - 3-d 

that generalizes the random walk behavior a = 1/2 of the d = 1 case. The physical meaning of the relation H24|l 
is the following [sj 33] : optimal configurations whose end-points are separated by a distance r typically merge 
at a distance I ^ r^/" and are then identical up to the origin, so that their difference in free energy then scales as 
r^l'^ ^ For large distance r ^ however, the two paths meet only the origin, and thus experience statistically 
independent disorders. This is why the random potential saturates for large separation [3lll33 |. 



2. Large-scale thermal fluctuations 

In the droplet scaling theory [3ll l32| | , there are large-scale thermal fluctuations that come from the rare nearly 
degenerate 'states'. In our description with favorables tubes, if we consider the excitations that involve a length I of 
the polymer : there exists a large number N ^ P (|15|l of other less favorable tubes of energies l'^ with transverse 
distance A ~ with respect to the best tube. As a consequence, with pro bability T/l'^, there exists an excited tube 
with a free-energy difference of order T with respect to the best tube jsj]. The power-law distribution P(A) of the 
end-point transverse scale A of thermal excitations computed via dynamical field theory [32| can be understood as 
follows 

P{A)A'-'dA ^ ^6(1 < L) (26) 

with the corresponding size I ^ A^/"^, so that 

P(A) ~ ^u; ^-^^ with6 = d+^ =(i + 2- i (27) 

where the cut-off w comes from the finite-size of the polymer I < L. With the values obtained before H12ll3f) . we thus 
obtain in dimension < d < 2, the power-law exponent 

2 + 2d-d'^ 



that generalizes the exponent b{d = 1) = 3/2. 
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Chaos exponent 



The sensitivity to small changes in the disorder distribution |26l l36l l37j or to small changes in temperature |31j 
have been interpreted at a scaling level as follows : if the small random perturbation e induces a change of order eL" 
for the free-energy of the previous state, there is an instability if a > w for lengths larger than £c(e) ~ e""^ with the 
so called chaos exponent c = ^^3^ • For a random bond perturbation [2^ |3^ |33| , and for a temperature change [U , 
the exponent of the perturbation is in both cases a = 1/2. For the temperature change, this comes from the behavior 
of entropy fluctuations A5' ~ i^/^, that was conjectured in [sj and numerically checked [3lll39| . With the values 
obtained before fl^. the chaos exponent thus reads in dimension < d < 2 

2(4 -d) ^^^^ 

that generalizes the well-known exponent c{d ~ 1) = 1/6. Finally, let us mention that more general perturbations 
with other exponents a have also been studied in [371 140| , and that the scaling form for the free-energy decorrelation 
has been recently computed on a Berker lattice |4H. 



IV. GENERALIZATION TO OTHER DIMENSIONS AND/OR OTHER TRANSVERSE 

CORRELATIONS 

In this section, we generalize our approach to the case of a D-dimcnsional manifold embedded in a space of total 
dimension D -\- d ~ dt in the solid-on-solid approximation : the manifold is described by a d-dimensional vector field 
r{x), where x is the D-dimensional vector of internal coordinates. We consider the standard Hamiltonian 



where V{x, r) is a Gaussian random potential with correlator 



(30) 



V{x,r)V{x',r') = S^{x ~ x')C{\r - r'\) (31) 
where the asymptotic behavior of the correlation in transverse directions is governed by some exponent a 

C{r) cx r" (32) 

r — ^00 

For instance, the case of local disorder characterized by delta correlations also in transverse directions corresponds 
via scaling to the case a = —d, whereas a random- field disorder corresponds to a = 1. 

The dimensional analysis of the elastic term L^~^R^ shows that the free behavior in the absence of disorder is 
characterized by the exponent 

2 — D 

Vfree{D) = — ^ (33) 

that generalizes the random walk exponent VfreeiD = 1) = 1/2 of the polymer case. In the following, we will only 
consider the cases < D < 2 where I'free > 0. 

A. Local Imry-Ma argument with the elastic term 

The generalization of the free energy © reads 

/_ ~ T^Ij + u_l'^/\''J^ (34) 



For < 0, the minimization of the free energy with respect to r_ yields 



2 ^ 4 - D 
^4-0^1^- ^[^Ij ^_ _ ('35') 

4 - a ^ ' 
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and the corresponding scaling for the free energy reads 

/_ ~ u^^r- with u;_ = 2i._ - (2 - D) = + - D)« (3g) 

4 — a 

This solution is consistent as long as > I'free i-G- since < I? < 2 

2D 

-^<«<4 (37) 

The condition au_ < shows that the physical interpretation is very different according to the sign of a. For 
all a < 0, as in the special case of uncorrelated disorder a = —d, the solution obtained corresponds to u_ > 0, i.e. 
to a unfavorable region, where the global free energy is positive /_ > 0. On the contrary, for a > 0, the solution 
corresponds to u_ < i.e. to a favorable region where the global free energy is negative /_ < 0. 



B. Local Imry-Ma argument with the confinement term 

The generalization of the Imry-Ma argument Q with the confinement entropy reads 

D 

U-T\J^] -u^l^J^l'^ (38) 

(A detailed analysis of the confinement entropy can be found in |42j). For auj^ < 0, the minimization with respect to 
r+ yields 



r+^u:*''Ul+ with u+ ^ , , (39) 

+ + + + 4D + q;(2-D) ^ ' 

The corresponding scaling for the total free energy reads 

2D D [2D + (2 - D)a] 



/. - with ^ D - . "4D"+ (2-Dir (^0) 



This solution is consistent as long as < < Vfree, i.e. 

2D 



D 



< a (41) 



The condition au+ < shows that the physical interpretation is again very different according to the sign of a. For 
a < 0, as in the special case of uncorrelated disorder a = —d, the solution obtained corresponds to u+ > 0, i.e. to a 
favorable region, where the global free energy is negative /+ < 0. On the contrary, for a > 0, the solution corresponds 
to u+ < i.e. to an unfavorable region where the global free energy is positive /+ > 0. 



C. Global optimization for a < 



For a < 0, as in the polymer case with uncorrelated disorder, the optimal strategy will be to find a favorable region 
L_^_ ^ L with confinement ~ L'^+. To find the best region available, it will be worth to afford a global elastic 
cost of L^~^Rq of the same order of magnitude of the energy gain ~ of the favorable region. The balance 
between these two terms yields that the global wandering will be 

UJ+ + 2- D D(4: - D) + (2 - D)a , , 

i?G L" w th V = = — ^ V^-^ — ^ — '— (42) 

^ 2 AD+{2-D)a ^ ' 

The global free energy gain will be 



(43) 
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The subleading transverse length scale characterizing the confinement reads 

D(2-D) , , 

i?+ - with v+ = —, '— 44) 

+ + 4D + a(2-D) ^ ' 

Finally, the use of (|16|) for the random variable u+ yields, via the change of variables H40|) to the following decay for 
the probability distribution of the free energy F ~ 

V(F) ~ -r-r — eryT^) w th rj+ = '— (45) 

^ ' F^-oo \F\ \L'^ J 2D ^ ' 

D. Global optimization for a > 

For a > 0, the previous picture completely changes, because the favorable regions are not the confined solutions 
(+), but the swollen solutions (— ). As a consequence, the global optimization coincide with a solution (— ) : the 
global exponents for the transverse direction and the free energy are thus directly given by 

1/ ^1^- = (46) 

4 - a ^ ' 

2D + {2-D)a 

u = a;_ - '— 47 

4 — a 

Finally, the use of (|16|l for the random variable u_ yields, via the change of variables H36|l to the following decay for 
the probability distribution of the free energy E ^ E^ 

ViE) ^ l-(\^\\-mT with, = i^ (48) 

We now briefly describe our results with their domain of validity for the interesting cases of a polymer 13 = 1 or of 
interfaces d = 1, before we compare with the results of other methods. 

E. Results for the polymer with decaying correlations with exponent —2 < a < 

For a polymer D = 1 with correlations described by the exponent a (|32() with —2 < a < 0, the results are the same 
as for the polymer with uncorrelated disorder with the replacement d —a, i.e. there is confinement with exponent 

= (49) 

with the corresponding exponents for the free energy and the wandering 

2 + a 

uj ^l-2v+ = -—— (50) 
4 + a ^ ' 

V = —— 51 

4 + a 

with the tail exponent rj^ ~ 2 + a/2. The domain of validity in a is limited by a — > —2, where the confinement 
exponent reaches the free random walk exponent — > 1/2, and by a — > where the optimization mechanism changes. 

F. Results for the polymer with long-range correlations with exponent < a < 4 

For a polymer D = 1 with correlations described by the exponent a H32|) with < a < 4, there is no confinement, and 
the exponents are directly given by the exponents (— ) coming from the direct dimensional analysis of the Hamiltonian 

(52) 



4 - a 



7 + fy 

UJ ^2v-\^--— (53) 
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and the tail exponent 77 = (4 — a)/2. The comparison with (|51|l shows that the exponents are continuous at a = 
with the values = 0) = 3/4 and Lj{a = 0) = 1/2. Special interesting cases discussed in P are a = 1, and a — 2. 
In particular for the latter case a = 2 that corresponds effectively to the Larkin model, there exists an exact solution 
in terms of replicas yielding the exact free-energy distribution [43j : our approach is in agreement with the exact 
results for the global exponents v = 3/2, uj = 2 and for the exponent 77 = 1 of the tail of the free-energy distribution 



G. Results for interfaces with random bond disorder 



For the special cases of interfaces of internal dimension < Z) < 2 in a space of total dimension D + 1 = dt, 
the random bond disorder corresponds to uncorrelated disorder in the d — 1 transverse direction, i.e. the effective 
exponent is a = —d = — 1 < 0. The constraint H41|l leads to the domain of validity 2/3 < D < 2, i.e. in terms of the 
total dimension dt of space 



The confinement exponent reads 



l<dt<S (54) 



. (55, 



Note that it vanishes in the limit — > 3, because the free exponent Vfree also vanishes in this limit of a two-dimensional 
surface D ^ 2 and is replaced by logarithms, the physical meaning being that the confinement is not a big constraint 
anymore ^3 • The corresponding exponents for the free-energy and wandering reads 

. ^id.~md,~5) 

5dt - 7 

_7..-8-d? ^^^^ 



with the tail exponent 77+ = (Sdt — 7)/{2{dt — 1)). 



H. Results for interfaces with random field disorder 



For the special cases of interfaces of dimension < Z? < 2 in a space of total dimension D + 1 = dt, the random field 
disorder corresponds to the exponent is a = +1 > 0. The constraint (|37|) does not modify the original constraints 
< < 2, i.e. in terms of the total dimension dt of space 

1< dt < 3 (58) 

In this case, there is no confinement, and the exponents are given by the exponents (— ) , i.e. they reduce to the well 
known exponents coming from the direct dimensional analysis of the Hamiltonian 

= (59) 
UJ = (60) 
with the left tail exponent rj = 3/2, that was also found in [23|. 



I. Agreement with Zhang replica scaling analysis 

The exponents presented in this Section are in agreement with Zhang replica-scaling analysis (21l] . It is thus useful to 
describe briefiy the 'dictionary' between the two methods. Whereas the traditional approach with replicas consists in 
considering the limit n — s- with the possibility of Replica Symmetry Breaking, the replica-scaling analysis proposed 
by Zhang 21] consists in analyzing the replicated problem within a symmetric treatment of the n replicas with n large 
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(for instance n(n — 1) is replaced by n^). The leading orders in (L, n) of the moments are then interpreted as coming 
from the tail of the probability distribution of the free-energy. In Zhang analysis, the case a < corresponds to an 
attraction between replicas and leads to a bound state, whereas the case a > corresponds to a repulsion between 
replicas with no bound state. So the presence of a bound state for replicas corresponds in our approach to the presence 
of a confinement tube. In conclusion, for those who prefer to think in real space than in replica space, our approach 
provides an equivalent self-contained description in real space; and for those who prefer to think in replica space, we 
hope that the translation in real space is also interesting, and can perhaps be useful in other disordered models. 

J. Comparison with other methods 

Within the field theoretical framework, the correctness of the direct dimensional Imry-Ma analysis for the random- 
field disorder a = 1 was supported by a functional RG analysis |T^, because the large distance behavior of the 
correlations are not renormalized, whereas for random-bond disorder a = —1, there are non-trivial renormalization of 
the exponents jl7j . So there should be a critical value ac below which the direct dimensional analysis does not give the 
correct exponents. It has been argued in various RG methods ilSj, 111 H H that (i) below the critical value ac, the 
wandering exponent i' sticks to the random bond value prb and would thus be 'superuniversal' (ii) the critical value 
ac is strictly negative. However, since the argument by Fisher (TtI that the long-range correlations of the disorder are 
not renormalized, naturally stops to apply at a = 0, it seems to us that from this point of view, the simplest scenario 
is that the true critical value is actually ac = 0, in agreement with the prediction of the replica-scaling analysis pll | 
and with our analysis with confinement. Moreover, we believe that our description involving a confinement mechanism 
suggests that a more appropriate field theoretical description in the case of transverse decaying correlations a < 
should perhaps include the presence of two different important transverse length scales, and should perhaps involve 
instantons calculations or other non-perturbative tools to describe the confinement. 

Finally, we should mention that the various alternative predictions 0, 143. |45| for the exponents in the domain 
a < have been tested via numerical simulations on kinetic roughening with various conclusions |46l . We believe 
that our approach, that explains the meaning of Zhang replica-scaling exponents in physical space strongly suggests 
that these exponents are indeed the correct ones. 

V. CONCLUSION 

In this article, we have analyzed via Imry-Ma energy/entropy arguments the strong disorder phases that exist in 
low dimensions at all non-zero temperatures for directed polymers and interfaces in random media. Within the field 
of disordered systems, the originality of random manifolds is that they have some freedom to 'choose' the disorder 
variables they see in a given sample, in contrast with spin systems for instance that cannot avoid any disorder variables. 
Our main result is that they can use this freedom to follow two qualitatively different strategies : (i) for disorder 
with decaying transverse correlations a < 0, the optimal strategy consists in being confined in a wandering tube, 
i.e. there are two different transverse scales that are important, namely the wandering scale and the confinement 
scale; (ii) for disorder with long-range transverse correlations a > 0, the optimal strategy is to be swollen, i.e. there 
is only one important transverse scale, and the exponents are given by a simple dimensional analysis of the original 
Hamiltonian. For the general case of a manifold of internal dimension D in a space of total dimension D + d in the 
presence of transverse correlations of exponent a, our results for exponents agree with Zhang replica-scaling analysis 
|2l|. Our method thus describes the same physics, but in real space instead of being in replica space, and this allows 
to make the link with the scaling droplet approach : the 'states' at finite temperatures of the scaling droplet theory 
|3ll l32j are interpreted in our approach as Imry-Ma favorable tubes, where the polymer remains confined, and it is 
this confinement in an Imry-Ma favourable tube that translates into a bound state for replicas in the replica-scaling 
analysis [2ll| . 

In conclusion, we hope that our approach via simple scaling arguments in real space gives a clear insight into the 
strong disorder phases in low dimensions, and a complementary point of view with respect to the other methods. 
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APPENDIX A: INITIAL DISORDER DISTRIBUTIONS WITH ALGEBRAIC TAILS 

For the lattice uncorrelated models, where each link has an energy V{i, x) drawn with some symmetric law P{V) = 
P{—V), the question of the universality with respect to the form of the initial probability distribution P{V) has 
attracted a lot of interest , because it appeared that the presence of long tails in the microscopic disorder 

was able to change the exponents [iM , Isof , as well as the morphology of the associated ultrametric tree structure 
of the locally optimal paths j^, even if the variance was finite /i > 2. 



1. The strategy of finding the best site energy at zero-temperature 

For the disorder (|A1|I , the heuristic argument 0, ^ that has been proposed to explain these dependences in /i 
consists in the balance between the maximal site energy Vmax = {LR'^Y^^^ drawn in the volume {LR^) and the global 
elastic energy B? /L, that yields the exponent 



^iV^alM = ^^d ^^^^ 



This exponent corresponds to the strategy of finding the maximal site energy and is usually considered as a lower 
bound 50] for the true exponent For instance in d = 1, this exponent is bigger than the value vcauss = 2/3 for 

/i > 5, which means that, at least for /i > 5, the true exponent is not the usual one vcauss = 2/3. There has 
been a large number of numerical studies on the true exponent either on the directed polymer problem or in 

corresponding growth models : it is usually believed P| that the exponent (|A2() is not the exact one 48, ,5fl| and in 
particular that the associated critical value yLtc((i = 1) = 5 in d = 1 is too low |52| and should be replaced by at least 
/ic(d = 1) ~ 7, even if it has been also argued that these numerical conclusions were due to cross-over effects and that 
the exponent (jA2|l was exact [s^ . 

From our point of view, it is not very clear why the exponent (|A2|) should describe the exponents at zero temperature 
of the lattice models that are measured in the numerical simulations. Indeed, for lattice models at zero temperature, 
there is no elastic energy, since a path is either allowed or forbidden. It is thus useful to reconsider in this context 
the strategy of finding the maximal site energy V„iax among a number of order L^'^'^ of independent sites that can 
be reached for a path of length L. The maximal value is thus of order Vmax ~ and since the distribution of 

its position is uniform, it reads in terms of the longitudinal coordinate I and the transverse distance r in continuum 
notations 

r'^-^e(r <1<L) 

PLil,r) L_ L (A3) 

As a consequence, the probability distribution of the transverse distance r of the best site Vmax takes the scaling form 

P^ir)^£dlp,il,r)^'-{^y-\l^L) (A4) 
i.e. the exponent corresponding to the strategy of finding the best site energy is 

-!v:::;(M,rf) = i (as) 

and any deviation from this extremal value = 1 on the lattice has to come from some cooperative effect of a large 
number of subleading best sites. In the following, we will not discuss the zero-temperature anymore, but the finite 
temperature case, and we will try to describe the cooperative effects by reconsidering our previous Imry-Ma arguments 
to see what changes are necessary in the presence of algebraic tails. So we first need to consider the random energy 
that can be gained in a tube. 



2. Distribution of tiie energy (7 of a path of length Z in a tube of radius r 



Let us consider the random energy 



C/ = ^F(z,x(z)) (A6) 
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of a path of length I in a tube of radius r in dimensions (1 + d). It is clear that its moments diverge for k > /j, as the 
moments of the initial law HA1|I . The tail of its probability distribution will thus present the same power-law decay 

where the scaling of the prefactor c(/, r) with respect to (/, r) can be estimated as follows. Let us introduce temporarily 
a large cut-off A in the original distribution IjAll) to regularize the diverging moments k > fi 

< V'' >A ^ A''-" (A8) 
The largest divergence in the cut-off A of the moment 

< C/fc >A^J2-J2< V{H,x{i,))...Viik,x{ik)) >A (A9) 

il ik 

is then given by the term of order < V'' >a, corresponding to coincident indices ii — i2 — ■■ = ik and coinciding 
transversal positions = ... = Xi^. If one assumes a uniform confinement in transverse directions within a tube of 
radius r, we thus obtain 

i.e. the prefactor c{l, r) of the tail (jA7|l scales as 

^(^'0-^ (All) 

On the other hand, the largest contribution in I of the moments (IA9|I with even k (odd moments vanish by symmetry) 
comes from the term corresponding to {k/2) pairs of coinciding indices, leading to 



k/2 



The natural rescaling appropriate for this term corresponds as it should to the rescaled variable u introduced before 
for the Gaussian case (0) 

U = u^I (A13) 
However, the rescaled variable u will now present the following algebraic tail I|A7IA11|I 

V{u) ^ (A14) 

The presence of a non-trivial radius r thus generalizes the special case r ~ 1 corresponding to the classical problem 
of the sum U of / independent variables : 

(i) For fj, > 2, the weight {lr'^y~^ of the tail of the rescaled variable u = U/^/l vanishes in the limit I oo, i.e. 
there is a generalized Central Limit Theorem describing the limit law of the rescaled variable u as soon as the variance 
is finite. However, since the polymer will try to find the best tube available, the presence of the algebraic tail can 
induce a change in global exponents as we will see. 

(ii) For fj, < 2, the weight of the tail in (|A14|I is a now positive power of I : the variable u is not appropriate 
anymore, and the new appropriate rescaled variable v is then defined by (|A7IA1H) 

U^vLT^R'^^ (A15) 
so that the limit probability distribution of v presents the tail 

Since the exponent of the transverse scale R in ljA15|) changes of sign at /i = 1, we will obtain a qualitative change at 
/i=l. 
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3. Local Imry-Ma arguments 

In this Section, we discuss the 'local' Imry-Ma arguments in the different cases /x>2, 1<^<2 and finally 
< /i < 1. 

a. Local Imry-Ma arguments for /j. > 2 

For /i > 2, the typical energy gained in a tube is still defined in terms of the rescaled variable u ljA13IA14|l . and 
thus the exponents l|10lll|l for typical favorable regions and (ti;_,t/_) I|7I8|I for typical unfavorable regions 

are unchanged. 

6. Local Imry-Ma arguments for 1 < /j. < 2 

For /i < 2, the new appropriate rescaled variable is v (|A15p . and thus the typical energy associated to a tube 
is different from the Gaussian case. The Imry-Ma arguments for favorable and unfavorable regions will thus yield 
different typical exponents. 

(i) Analysis of the exponents in typical favorable regions 

The free energy of a typical region (?_(_, r_|_) with the confinement term then becomes 

(A17) 

For a favorable region > 0, the minimization of /+ with respect to r_|_ yields 
The corresponding free-energy reads 

u ^ -.r*-,-"> W.U, - ^fii^ (AW) 

For < d < 2, this analysis is valid for 1 < /x < 2 : the confinement exponent varies between — s- 1) — s- and 
its usual value 2) 1/(4 — d), and the free-energy exponent varies between — > 1) ^ 1 and its usual 

value — !• 2) ^ (2 — d)/(4 — c?). In particular, in dimension d ~ I, the exponents read 

=^ (A20) 

u;+(M,d=l) =^ (A21) 
+ 1 

For d > 2, these exponents are expected to be valid as long as the confinement exponent is positive and 

smaller than the free value 1/2, i.e. for ji < 1 -\-2/d For instance, in transverse dimension d ~ 2, the validity domain 
is 1 < /X < 2, whereas inc? = 3itisl</Lt<5/3. Equivalently, for fixed /i in the interval 1 < /Lt < 2, the exponents 
are valid for d < 2/ {p. — 1). 

(ii) Analysis of the exponents in typical unfavorable regions 

The free energy of a typical region (Z_,r_) with the elastic term © then becomes 

f_ ^T— + v^Vlry' ^ (A22) 

For a unfavorable region u_ > 0, the minimization of /_ with respect to r_ yields 

with ,y_(^) = — ii±i — - (A23) 
{2 + d)^-d ^ ' 



The corresponding free-energy reads 



-(M) ,., 2 + d-dii 



/_ ^ .™-r_- with . ^1-^ (A24) 
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For < c? < 2, this analysis is valid for 1 < ^ < 2 : the roughness exponent varies between — > 1) — > 1 and its 
usual value V-{p. ^ 2) — > 3/(4 + d), and the free-energy exponent varies between ^-(/i — *■ 1) — > 1 and its usual value 
ti>_(/i ^ 2) ^ (2 — (i)/(4 + d). In particular, in dimension d—1, the exponents read 

v.{^Ji,d=l) (A25) 

c._(/i,d=l) =1^3^ (A26) 



c. Local Imry-Ma argument for < /j. < 1 



For < < 1, the exponent of the transverse scale R in (jA15|) is positive. As a consequence, exactly as in the 
case of long-range transverse correlation, the favorable regions do not correspond to confined solutions, but to swollen 
solutions. The free-energy of a region (Z_,r_) with the elastic term © has now to be replaced by 

/_ -T— (A27) 

For a favorable region u_ < 0, the minimization with respect to r„ yields 



{2 + d)^-d ^ ' 

and the corresponding free-energy reads 

In particular, in the limit /i 1, the free-energy exponent — > 1^) = 1 is in continuity with the 0J+{^i ^ 1^) = 1 



4. Global optimization for the full polymer 



In the presence of long tails, the global optimization for the full polymer consists in looking for the confinement 
scale Rs ~ L'^s(m) ^j^cJ the global wandering scale Rq L^'^>^^ that minimize the global free energy containing three 
contributions 

Fl{Rg, Rs) = + T A „ U,rr.ax{L, Rs;N = §) (A30) 

s s 

The first term Rq/L represents the global elastic cost for the diffusive tube. The second term represents the 
confinement entropy cost. The third term Umax represents the best energy among N independent variables U, where 
U represents the random energy associated to a tube (i, Rs)- Using the form 1A7IA1H) of the tail of its probability 
distribution, we thus obtain the following scaling 

UmaAL, Rs; N^^)^ (c(L, Rs)N)^ ^ i?^'^ (Li?^) - (A31) 

Rg 

This result shows why the introduction of two transverse scales (Rs, Rg) is actually always better to gain energy than 
the swollen solution that involves only one scale Rq, whose free-energy reads 

Fr°"^"(i?G) = - UtypiL, Rg) (A32) 

where the typical value scales as Utyp{L, Rg) ^ Rq'^ (Li?^) for < 2 HA15() . The comparison with HA31|I is thus 
immediate : both have the same factor {LR'q) ^ presenting the scale of the best sites in the volume {LRq), but the 
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factor of the density in the transverse directions of the tube is much better in a confined tube Rg than for the swollen 
solution Rq"^- 

The optimization of ljA30|) with respect to {Rs,Rg) yields the following exponents 

-'^) (-^) 

These exponents are valid in transverse dimension < d < 2 in the domain 

^s(rf) = + 1 < M < Mc(rf) (A35) 

where the inferior value fj,s{d) = d+l corresponds to the point where the confinement exponent vanishes vsid+l) = 0, 
whereas the superior value corresponds the critical value 

where the Gaussian exponents are recovered. This critical value varies between udd — > 0) ^ 2 and fic{d ^ 2) — > 2, 
and is maximal for d = 1 with the value udd = 1) = 5/2. For instance, in dimension d = 1, the exponents (|A34|I 

ly{^i,d=l) =^ (A37) 

Mt^,d^l) =^ (A38) 

are thus expected to be valid for 2 < fi < fic{d = 1) = 5/2. 

Let us now discuss what happens for ii < fis{d) — d+ I. Below this value, the confinement exponent would become 
negative, which is unphysical, and thus the confinement exponent will stick to its minimal value 

lysifi < ^^sid) = d+l) = (A39) 

corresponding to an extreme confinement. In this regime, the confinement entropy of order L is subleading with 
respect to the elastic entropy and the disorder energy in ljA30|) . The balance between these two terms then coincides 
with the estimation (|A2|I coming from the single best site strategy 



^(^ < < ^,sid) =d+l) = 4Tnt)(M) = (A40) 



Note that this wandering exponent is then always greater than one. As a consequence, on a lattice, the wandering 
exponent will stick to its maximal value v = I 

^latUce^^ < /i,(d) = d + 1) = 1 (A41) 

5. Discussion 

In this Appendix, we have shown how the presence of algebraic tails in the initial distribution could be taken into 
account in our approach to yield a change in the global exponents even in the cases where the variance is finite fi > 2. 
However, even if our results are on this point qualitatively correct, the status of the quantitative results is not clear. 
Indeed, in dimension d = 1, we have obtain as fi varies for the wandering exponent at non-zero temperature 

z/(d= 1,2 < ^ < 5/2) =^—- (A42) 

v(d = 1, 1/2 < < 2) = Ji±l. > 1 i.e. < 2) = 1 (A43) 

2fi — 1 

where the change that takes place at ^ = 2 with the corresponding value v{ix = 2) = 1, between a phase with a positive 
confinement exponent Vs{^^ > 2,d = 1) > and a phase with a vanishing confinement exponent Vsip, < 2,d = 1) = Q. 
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The results (|A43|) for > 2 are far from the exponents measured at T = in numerical simulations : in particular, the 
critical value we have obtained fic{d = 1) = 5/2 is very far from the lattice numerical estimation udd = 1, T = 0) ~ 7 
0, Is there a problem in this case between the zero-temperature best path and our approach at non-zero 

temperature based on the confinement entropy? Or is our Imry-Ma approach too simple to describe correctly the 
algebraic tails ? For instance, it may be that the smaller scales have to be taken into account, in contrast with the 
Gaussian case where the optimization on the biggest scale fixes the exponents. 
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